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Abstract. Various origins of linear and nonlinear Schrodinger equations are discussed in 
connection with diffusion, hydrodynamics, and fractal structure. The treatment is mainly 
expository, emphasizing the quantum potential, with a few new observations. 
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1. INTRODUCTION 

Perhaps no subject has been the focus of as much mystery as "classical" quantum me- 
chanics (QM) even though the standard Hilbert space framework provides an eminently 
satisfactory vehicle for determining accurate conclusions in many situations. This and 
other classical viewpoints provide also seven decimal place accuracy in QED for example. 
So why all the fuss? The erection of the Hilbert space edifice and the subsequent de- 
velopment of operator algebras (extending now into noncommutative (NC) geometry) has 
an air of magic. It works but exactly why it works and what it really represents remain 
shrouded in ambiguity. Also geometrical connections of QM and classical mechanics (CM) 
are still a source of new work and a modern paradigm focuses on the emergence of CM 
from QM (or below). Below could mean here a micro structure of space time (quantum 
foam, Cantorian spacetime, etc.). In addition there are beautiful stochastic theories for 
diffusion and QM. In terms of background information in book form we mention here e.g. 

laiiniiiaEiiEaEHiEiiiMiEaizziizniEaiHaEsiEEi^ (the 
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lecture notes ^1 El ^1 ^1 E| in a more polished and organized form should also eventu- 
ally become part of a book in preparation). The present paper focuses on certain aspects 
of the Schrodinger equation (SE) involving the wave function form ■0 = Rexp{iS/h), hydro- 
dynamical versions, diffusion processes, quantum potentials, and fractal methods. The aim 
is to envision "structure" , both mathematical and physical, and we avoid detailed technical 
discussion of mathematical fine points (cf. [HI |2H1 El IM1 EH Ell CSl USH] for various 
delicate matters). Rather than looking at such matters as Markov processes with jumps for 
example we prefer to seek "meaning" for the Schrodinger equation via microstructure and 
fractals in connection with diffusion processes and kinetic theory. 



2. BACKGROUND FOR THE SCHRODINGER EQUATION 

First consider the SE in the form (Al) — (h"^ /2m)'tp" + Vtp = ihijjt so that for ij) 
Rexp{iS/h) one obtains 



2 u" 



(2.1) 



Im Zmn 



1 



0; dt{R') + -iR'Sy = 



m 



where S' ~ dS/dX. Writing P = R^ (probability density ~ 
(quantum potential) this becomes 

(2.2) St + ^ + Q + U = 0; Pt + -{PS')' = 

2m m 



andQ = -{h^/2m){R" /R) 



and this has some hydrodynamical interpretations in the spirit of Madelung. Indeed going 
to for example we take p = S' with p = mq for q a velocity (or "collective" velocity - 
unspecified). Then 1)2. 2(1 can be written as {p = mP is an unspecified mass density) 
(2.3) 

St + f + Q + V = 0; Pt + -{Pp)' = 0;p = S'; P = R^- Q = = 

2m m 2m R 2m ^ 



Note here 
(2.4) 



_ 1 



VP 4 

Now from S' = p = mq = mv one has 



p 



(2.5) 



P, + {Pqy = 0^p, + ipqy = 0; St+^ + V-^^=0 

2m 2m p 



Differentiating the second equation in X yields {d ~ d/dX, v = q) 



(2.6) 



mvt + mvv' + dV 



2m \ J 

Consequently, multiplying by p = mv and p respectively in (|2.5() and ()2.6|) . we obtain 



(2.7) mpvt + mpvv' + pdV pd ( — 

2m V \/P 



0; mvpt + mv{p'v + pv') = 
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Then adding in (|2.7|) we get 

(2.8) dtipv) + dipv') + ^dv - -^pd ( ) = 

This is similar to an equation in |1J (called an "Euler" equation) and it definitely has a 
hydrodynamic flavor (cf. also '60'). 

Now go to |124j and write 1)2. 6(1 in the form {mv = p = S') 

3v 1 

(2.9) — + {vV)v = V{V + Q); vt + vv' = -{l/m)d{v + Q) 

ot m 

The higher dimensional form is not considered here but matters are similar there. This 
equation (and (|2.8() ) is incomplete as a hydrodynamical equation as a consequence of a 
missing term —p^^Vp where p is the pressure (cf. [HH)- Hence one "completes" the equation 
in the form 

(2.10) m + {v ■ V)v^ = -V{V + Q)- VF; mvt + mvv' = -d{V + Q) - F' 

where (A2) VF = (l/ii^)Vp (or F' = (l/i2^)p'). By the derivations above this would then 
correspond to an extended SE of the form 

(2.11) = -—Ai^ + ViP + FiP 

ot 2m 

provided one can determine F in terms of the wave function ip. One notes that it a nec- 
essary condition here involves curlgrad{F) = or (A3) curl{R^^\/p) = which enables 
one to take e.g. (A4) p = —bB? = —h\'ijj\'^ . For one dimension one writes (A5) F' = 
-5(l/i?2)a|^|2 = -{2bR'/R) ^ F = -2blog{R) = -blog{\^\'^). Consequently one has a 
corresponding SE 

(2.12) = -^r + - fc(%iV'p)V' 

ot 2m 

This equation has a number of nice features discussed in jl24j (but serious drawbacks as 
indicated in 23, - cf. also [SZl IHl 1121 113 EBl ITTTTl ITTM Tm\ \. For example (A6) ^ = 
f5G{x — vt)exp{ikx — iiot) is a solution of (|2.12j) with F = and for v = hk/m one 
gets (A7) ijj = cexp[—{B / A){x — vt + d)'^]exp{ikx — itot) where B = Amb/ti? . Normal- 
ization = 1 is possible with (A8) l^/^p = 5m{Cj = sj maj'nexp^—am^) where 
a = 2b/fi^ , d = 0, and ^ = x — vt. For m — > oo we see that 6m becomes a Dirac delta and 
this means that motion of a particle with big mass is strongly localized. This is impossible 
for ordinary QM since exp{ikx — iut) cannot be localized as m — > oo. Such behavior helps 
to explain the so-called collapse of the wave function and since superposition does not hold 
Schrodinger's cat is either dead or alive. Further v = kh/m is equivalent to the deBroglie 
relation A = h/p since A = {2iT/k) = 2iT{h/mv) = 27r(/i/27r)(l/p). 

REMARK 2.1. We go now to [ZOl and the linear SE in the form (A9) i{d'4)/dt) = 

— (l/2m)A'0 + C/(r)^; such a situation leads to the Ehrenfest equations which have the form 
(AlO) < V >= {d/dt) <f> and < f >= / (i^x|^/^(r, t)pr and (All) m{d/dt) < v >= F{t) 
with Fit) = — J d'^x\ip{f,t)\'^'VU{f). Thus the quantum expectation values of position and 
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velocity of a suitable quantum system obey the classical equations of motion and the ampli- 
tude squared is a natural probability weight. The result tells us that besides the statistical 
fluctuations quantum systems posess an extra source of indeterminacy, regulated in a very 
definite manner by the complex wave function. The Ehrenfest theorem can be extended 
to many point particle systems and in [TOj one singles out the kind of nonlinearities that 
violate the Ehrenfest theorem. A theorem is proved that connects Galilean invariance, and 
the existence of a Lagrangian whose Euler-Lagrange equation is the SE, to the fulfillment 
of the Ehrenfest theorem. ■ 

REMARK 2.2. There are many problems with the quantum mechanical theory of 
derived nonlinear SE (NLSE) but many examples of realistic NLSE arise in the study of su- 
perconductivity, Bose-Einstein condensates, stochastic models of quantum fluids, etc. and 
the subject demands further study. We make no attempt to survey this here but will give 
an interesting example later from [2Sj related to fractal structures where a number of the 
difficulties are resolved. For further information on NLSE, in addition to the references 
above, we refer tolZlEHlElEElIZnilZlllZSl 11221 EHOl CSll 11101 Ell for some typical situa- 
tions (the list is not at all complete and we apologize for omissions). Let us mention a few 
cases. 

• The program of (THI introduces a Schrodinger Lagrangian for a free particle includ- 
ing self-interactions of any nonlinear nature but no explicit dependence on the space 
of time coordinates. The corresponding action is then invariant under spatial coor- 
dinate transformations and by Noether's theorem there arises a conserved current 
and the physical law of conservation of linear momentum. The Lagrangian is also re- 
quired to be a real scalar depending on the phase of the wave function only through 
its derivatives. Phase transformations will then induce the law of conservation of 
probability identified as the modulus squared of the wave function. Galilean invari- 
ance of the Lagrangian then determines a connection betwee the probability current 
and the linear momentum which insures the validity of the Ehrenfest theorem. 

• We turn next to ^ for a statistical origin for QM (cf. also ;llti38l[7TT[73l[lTn[TT8l 
I133j ). The idea is to build a program in which the microscopic motion, underlying 
QM, is described by a rigorous dynamics different from Brownian motion (thus 
avoiding unnecessary assumptions about the Brownian nature of the underlying 
dynamics). The Madelung approach gives rise to fluid dynamical type equations 
with a quantum potential, the latter being capable of interpretation in terms of a 
stress tensor of a quantum fluid. Thus one shows in ^2] that the quantum state 
corresponds to a subquantum statistical ensemble whose time evolution is governed 
by classical kinetics in the phase space. The equations take the form 



(2.13) pt + 9,(H = 0; dtiiJipui) + 5j(p0,,) + pd^y = 0; dt{pE) + d^{pS) - pdtV = 



(2.14) 
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(2.15) 



for two scalar fields S determining a quantum fluid. These can be rewritten as 

dt 



1 d'^S }_d^dS__^ 
fj, dx"^ fj, dx dx 



(2.16) 



4/i dx"^ 

where ^ = log{p) and for Vt 
h = Nr] one arrives at a SE 

ih 



= (e/2) + {il^)S 



+ V = Q 



dx 

log^ with m 



Nfi, V = NV, and 



dt 2m dx"^ 

Further one can write ^I* = p^^'^exp{i&/h) with & = NS and here N = J \^\'^(Px. 
The analysis is very interesting. ■ 

REMARK 2.3. Now in [44j one is obliged to use the form ip = Rexp(iS/h) to make 
sense out of the constructions (this is no problem with suitable provisos, e.g. that S is 
not constant - cf. [HI CH EZl |1H| ) . Thus note from (A12) tp'/ip = (R'/R) +i{S'/h) with 
^{ip'/il:) = {l/m)S' ~ p/m (see also ^(TV^ below). Also note (A13) J = {h/m)^^!*^;' and 
p = R? = IV'P represent a current and a density respectively. Then using p = mv = mq one 
can write (A14) v = {h/m)Q{'ip' /ip) and J = {h/m)Q\^\^{^p*^' /{ipl^) = {h/m)Q{pv). Then 
look at the SE in the form ihipt = —{h?/2m)ip" + Vip with ipt = {Rt + iStR/h)exp{iS/h) 
andip^^ = [{R' + {iS'R/H)exp{iS/h)]' = [R" + {2iS' R' /h) + {iS" R/h) + {iS' /hfR]exp{iS/H) 
which means 



(2.17) 



2m 



R" 



2iS'r' iS"R 



h 



+ VR = ih 



Rt + 



iStR 



h^R" 



+ 1^ = 



This can also be written as 

(2.18) dtp + -d{pp) = 0; St + ^ + Q + V = Q 

m 2m 

where Q = —fi?R"/2mR. Now we sketch the philosophy of |441 145j in part. Most of such 
aspects are omitted and we try to isolate the essential mathematical features. First one 
emphasizes configurations based on coordinates whose motion is choreographed by the SE 
according to the rule (1-D only here) 



(2.19) 



V = — ^ , , ,^ 



m 



where (A15) ihipt = —{h'^ /2m)ip" + Vip. The argument for 1)2. 19(1 is based on obtain- 
ing the simplest Galilean and time reversal invariant form for velocity, transforming cor- 
rectly under velocity boosts. This leads directly to ()2.19|) (~ (A14)) so that Bohmian 
mechanics (BM) is governed by ()2.19|) and (A15). It's a fairly convincing argument 
and no recourse to Floydian time seems possible (cf. ^2 HHl EDI 1^)- Note however 
that if 5 = c then q = v = {h/m)Q{R' /R) = while p = 5' = so perhaps this 
formulation avoids the 5 = problems indicated in jlll 1481 I5UI I51j . One notes also 
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that BM depends only on the Riemannian structure g = (gij) = {rriidij) in the form 
(A16) q = hQ{gradtp/ilj); ihipt = — (^^/2)A^ + Vip. What makes the constant h/m 
in (|2.19j) important here is that with this value the probability density IV'P on config- 
uration space is equivariant. This means that via the evolution of probability densities 
pt + div{vp) = (as in (|2.18)) with v ~ p/m) the density p = I'i/'p is stationary relative 
to ip, i.e. p{t) retains the form \ip{q,t)\'^ . One calls p = IV'P the quantum equilibrium 
density (QED) and says that a system is in quantum equilibrium when its coordinates are 
randomly distributed according to the QED. The quantum equilibrium hypothesis (QHP) 
is the assertion that when a system has wave function ip the distribution p of its coordinates 
satisfies p = |V'P- ■ 

REMARK 2.4. We extract here from [^IHHinSl (cf- also the references there for back- 
ground and IHH] for some information geometry). There are a number of interesting 
results connecting uncertainty, Fisher information, and QM and we make no attempt to 
survey the matter. Thus first recall that the classical Fisher information associated with 
translations of a 1-D observable X with probability density P{x) is 

(2.20) j dxPix){[log{Pix)]'f > 

One has a well known Cramer-Rao inequality (A17) Var{X) > F^^ where Var{X) ~ 

— 1/2 

variance of X. A Fisher length for X is defined via (A18) 6X = and this quantifies 

the length scale over which p{x) (or better log{p{x))) varies appreciably. Then the root 
mean square deviation /S.X satisfies (A19) AX > 5X. Let now P be the momentum 
observable conjugate to X, and P^i a classical momentum observable corresponding to the 
state tp given via (A20) Pci{x) = {h/2i)[{^' /tp) - (tp'/ip)] (cf. ^J^). One has the identity 
(A21) < p >^=< Pel >^ following from (A20) with integration by parts. Now define the 
nonclassical momentum by pnc = P — Pel and one shows then (A21) AXAp > 5X^p > 
5X^Pnc = h/2. Now go to now where two proofs are given for the derivation of the SE 
from the exact uncertainty principle (as in (A21)). Thus consider a classical ensemble of 
n-dimensional particles of mass m moving under a potential V. The motion can be described 
via the HJ and continuity equations 

(2^21) | + J_|v,p + K = 0;f + V.[f^1=0 

ot 2m at \_ m 

for the momentum potential s and the position probability density P (note that we have 
interchanged p and P from [S^] - note also there is no quantum potential and this will be 
supplied by the information term). These equations follow from the variational principle 
6L = Q with Lagrangian 



(2.22) L = dtdJ'xP 



ds ^ 1 1^ ,2 _|_ 
dt 2m 

It is now assumed that the classical Lagrangian must be modified due to the existence of 
random momentum fluctuations. The nature of such fluctuations is immaterial for (cf. [HI] 
for discussion) and one can assume that the momentum associated with position x is given 
by (A22) p = Vs + A'^ where the fluctuation term N vanishes on average at each point x. 
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Thus s changes to being an average momentum potential. It follows that the average kinetic 
energy < |Vsp > /2m appearing in (|2.22|) should be replaced by < |Vs + A^p > /2m giving 
rise to 

(2.23) L' = L + (2m)-i j dt < N ■ N >= L + (2m)"^ j dt{ANf 

where AA^ =< N ■ N >^/^ is a measure of the strength of the fluctuations. The addi- 
tional term is specified uniquely, up to a multiplicative constant, by the following three 
assumptions 

(1) Action principle: L' is a scalar Lagrangian with respect to the fields P and s where 
the principle 5L' = yields causal equations of motion. Thus (A23) (AA^)^ = 
/ d'^xpf{P,VP,dP/dt,s,'Vs,ds/dt,x,t) for some scalar function /. 

(2) Additivity: If the system comprises two independent noninteracting subsystems with 
P = P1P2 then the Lagrangian decomposes into additive subsystem contributions; 
thus (A24) / = /i + /2 for P = P1P2. 

(3) Exact uncertainty: The strength of the momentum fluctuation at any given time 
is determined by and scales inversely with the uncertainty in position at that time. 
Thus (A25) AA^ — > kAN for x x/k. Moreover since position uncertainty is 
entirely characterized by the probability density P at any given time the function / 
cannot depend on s, nor explicitly on t, nor on dP/dt. 

The following theorem is then asserted (see |131 162j for the proofs) . 

THEOREM 2.1. The above 3 assumptions imply (A26) (AA^)^ = cj d'^x P\Vlog{P)\^ 
where c is a positive universal constant. 

COROLLARY 2.1. It follows from ()2.23|) that the equations of motion for p and s corre- 
sponding to the principle 6L' = are 

(2.24) ,,^._|lv=^ + V« 

where h = 2y/c and ^ = VPexp{is/h). ■ 

REMARK 2.5. We sketch here for simplicity and clarity another derivation of the SE 
along similar ideas following jl32j . Let P{y^) be a probability density and P{y^ + Ay*) be 
the density resulting from a small change in the y*. Calculate the cross entropy via 

(2.25) J(P(y* + Ay') : P{y')) = j P{y' + Ay^og^^^^^^^^d'^y ^ 
If 1 dP{y')dP{y' 



-ry 



Ay* A/ = /,fcAy*A/ 



2 J P{y^) dy' dy^) 
The Ijk are the elements of the Fisher information matrix. The most general expression 
has the form 



(ooa\ T ^ / 1 dP{x'\e')dP{x'm 

where P{x'^\9'^) is a probability distribution depending on parameters 9' in addition to the 
x\ For (A27) P(x*|6'*) = P(x* + 6^) one recovers (IT^ (straightforward - cf. [HI]). If 
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P is defined over an n-dimensional manifold with positive inverse metric g^^ one obtains a 
natural definition of the information associated with P via 

(2.27) 



2 J Pdy^dy^ 

Now in the HJ formulation of classical mechanics the equation of motion takes the form 
(2.28) BS I BS8S 

where g'^'^ = diag{l/m, • • • ,1/m). The velocity field n'^ is given by (A28) n'^ = g'^'^{dS/dx'^). 
When the exact coordinates are unknown one can describe the system by means of a proba- 
bility density P{t, with (A29) / Pd^'x = 1 and (A30) {dP/dt)+{d/dxf){Pgf^{dS/dx'') -- 
0. These equations completely describe the motion and can be derived from the Lagrangian 



^ dt 2 ^x^^ dx 

using fixed endpoint variation in S and P. Quantization is obtained by adding a term 
proportional to the information I defined in H2.27() . This leads to 

f r.^^S 1 dS dS X dP dP] ^1 , ,„ 

Fixed endpoint variation in S leads again to (A30) while variation in P leads to 

+ V = 



(2.31) ^ + -g'^'^ 

^ ^ dt 2^ 



dS dS ^ f I dP dP 2 d'^P 
+ A 



dx^^ dx'^ \ P2 dx^^ dx^ P dx^^dx^ 



These equations are equivalent to the Schrodinger equation if (A31) ip = y Pexp{iS /h) 
with A = {2hf (cf. Section 6). ■ 

REMARK 2.6. The SE gives to a probability distribution p = (with suitable 
normalization) and to this one can associate an information entropy S{t) (actually config- 
uration information entropy) (A32) S = — J plog{p)d^x which is typically not a conserved 
quantity (S is an unfortunate notation here but we retain it momentarily since no confusion 
should arise) . The rate of change in time of S can be readily found by using the continuity 
equation (A33) dtp = —V • {vp) where v is & current velocity field Note here (cf. also [l26j 

d S f f 
(2.32) _ = _ y + log{p))dx = y (1 + log{p))d{vp) 

Note that a formal substitution of u = —u in (A33) implies the standard free Browian 
motion outcome (A34) dS/dt = D ■ j[{Vpf/p)d^x = D ■ Tr^ > - use (A35) u = 
DVlog{p) with D = h/2m) and (lO^ with / {l + log{p))d{vp) = -J vpdlog{p) = - / t-p' ~ 
jiip'Y'/p) modulo constants involving D etc. Recall here mfF ~ —{2/D'^) J pQdx = 
J (ix[(V/9)2/p] is a functional form of Fisher information. A high rate of information entropy 
production corresponds to a rapid spreading (flattening down) of the probablity density. 
This delocalization feature is concomitant with the decay in time property quantifying 
the time rate at which the far from equilibrium system approaches its stationary state of 
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equilibrium (A36) d/dtTr^ < 0. ■ 

REMARK 2.7. Now going back to the quantum context one admits general forms 
of the current velocity v. For example consider a gradient field v = b — u where the 
so-called forward drift b{x, t) of the stochastic process depends on a particular diffusion 
model. Then one can rewrite the continuity equation as a standard Fokker-Plank equation 
(A37) dtp = DAp — V • (bp). Boundary restrictions requiring p, vp, and bp to vanish at 
spatial infinities or at boundaries yield the general entropy balance equation 

dS 



d^x = -D-j^ = I p{v ■ u)d^x =< V ■u> 



P 

The first term in the first equation is not positive definite and can be interpreted as an 
entropy flux while the second term refers to the entropy production proper. The fiux term 
represents the mean value of the drift fleld divergence V • b which by itself is a local measure 
of the flux incoming to or outgoing from an infinitesimal surrounding of x at time t. If 
locally (V • b){x,t) > on an infinitesimal time scale we would encounter a local entropy 
increase in the system (increasing disorder) while in case (V • 6)(x, t) < one thinks of local 
entropy loss or restoration or order. Only in the situation < V • 6 >= is there no entropy 
production. Quantum dynamics permits more complicated behavior. One looks first for 
a general criterion under which the information entropy (A32) is a conserved quantity. 
Consider (|2.8j) and invoke the diffusion current to write (recall u = D(Vp)/p) 

(2.34) ^lt=~ j t'^^'^'^'^^)] • \p'"\DVp)]d\ 

Then by means of the Schwarz inequality one has (A38) D\dS/dt\ << >^/'^< v? y-^l"^ so 
a necessary (but insufficient) condition for dS/dt 7^ is that both < > and < > are 
nonvanishing. On the other hand a sufficient condition for dS/dt = is that either one of 
these terms vanishes. Indeed in view of (A39) < V? >= J[{Vpf/p]d^x the vanishing 
information entropy production implies dS/dt = 0; the vanishing diffusion current does the 
same job. ■ 

REMARK 2.8. We develop a little more perspective now (following - first paper). 
Recall Q written out as 



(2.35) Q = 20"^^^— = D'^ 

p^l^ 



1 9 



p 2p2 

where u = DVlog{p) is called an osmotic velocity field. The standard Brownian motion 
involves v = —u, known as the diffusion current velocity and (up to a dimensional factor) 
is identified with the thermodynamic force of diffusion which drives the irreversible process 
of matter exchange at the macroscopic level. On the other hand, even while the thermody- 
namic force is a concept of purely statistical origin associated with a collection of particles, 
in contrast to microscopic forces which have a direct impact on individual particles them- 
selves, it is well known that this force manifests itself as a Newtonian type entry in local 
conservation laws describing the momentum balance; in fact it pertains to the average (lo- 
cal average) momentum taken over by the particle cloud, a statistical ensemble property 
quantified in terms of the probability distribution at hand. It is precisely the (negative) 
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gradient of the above potential Q in (|2.35|) which plays the Newtonian force role in the 
momentum balance equations. The second analytical expression of interest here involves 

(2.36) - J Qpdx = (1/2) J u^pdx = (1/2)1)2 . ^ j ^^^dx 

where Fx is the Fisher information, encoded in the probability density p which quantifies 
its gradient content (sharpness plus localization/disorder) (note — / Qp = — f[{l/2)u'^ p + 
Dpu'] = -J{l/2Wp + fDup' = -{l/2)j-D^{p'/pfp + D^^p'{p'/p) = {D^ /2) f {p'f / p = 
(1/2) j u^p)- On the other hand the local entropy production inside the system sustaining 
an irreversible process of diffusion is given via 

(2.37) f = D.p^dx = D.Fx>0 

This stands for an entropy production rate when the Fick law induced diffusion current 
(standard Brownian motion case) j = —DVp, obeying dtp + Vj = 0, enters the scene. Here 
S = — J plog{p)dx plays the role of (time dependent) information entropy in the nonequi- 
librium statistical mechanics framework for the thermodynamics of irreversible processes. 
It is clear that a high rate of entropy increase coresponds to a rapid spreading (flattening) of 
the probability density. This explicitly depends on the sharpness of density gradients. The 
potential type Q(x,t), the Fisher information Fx, the nonequilibrium measure of entropy 
production dS/dt, and the information entropy S{t) are thus mutually entangled quantities, 
each being exclusively determined in terms of p and its derivatives. 

In the standard statistical mechanics setting the Euler equation gives a prototypical mo- 
mentum balance equation in the (local) mean 

F VP 

(2.38) {dt + vV)v = 

m p 

where F = —VF represents normal Newtonian force and P is a pressure term. Q appears 
in the hydro dynamical formalism of QM via 

1 1 A 0^/2 

(2.39) {dt + v V)v = -F-VQ = -F + —-V-P— 

Another spectacular example pertains to the standard free Brownian motion in the strong 
friction regime (Smoluchowski diffusion), namely 

(2.40) {dt + v- V)v = -2D2v-^ = -VQ 

p^i'^ 

where v = —D{SJp/p) (formally D = h/2rn). ■ 

REMARK 2.9. The papers in [5^] contain very interesting derivations of Schrodinger 
equations via diffusion ideas a la Nelson, Markov wave equations, and suitable "applied" 
forces (e.g. radiative reactive forces). ■ 
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3. DIFFUSION AND FRACTALS 

We go now to Nagasawa jHHl lEHl IHOl 1^ to see how diffusion and the SE are really 
connected (cf. also UHl iSl EH ESI EH HIH ESI IHZl 11201 HH] for related material, 
some of which is discussed later in detail); for now we simply sketch some formulas for a 
simple Euclidean metric where (Bl) A = ^(9/9x*)^. Then ^(t,x) = exp[R{t,x) + iS{t,x)] 
satisfies a SE (B2) idt^J + (1/2) A^p + ia{t, x) • VV' - V{t, x)^; = [h and m omitted) if and 
only if 

(3.1) v = -^+^-AR+ ^-{VRf - ^{VSf - a • V5; 

an 1 

= — + -AS + (VS) ■ (VR) + a-VR 

in the region (B3) D = {{s,x) : ■ip{s,x) ^ 0}. Solutions are often referred to as weak or 
distributional but we do not belabor this point. From |88L l9Uj there results 

THEOREM 3.1. Let 'tpit,x) = exp[R{t,x) + iS{t,x)] be a solution of the SE (B2); then 
(B4) (p{t, x) = exp[R{t, x) + S{t, x)] and (j) = exp[R{t, x) — S{t, x)] are solutions of 

(3.2) ^ + ^Ac/} + a{t, x)-V(p + c(t, x, (/))0 = 0; 

+ - a(t, x)-V^ + c(t, X, (t))4) = 
where the creation and annihilation term c{t,x, 0) is given via 

QQ 

(3.3) c(t, X, (j)) = -V{t, x) - '^-Q^it, x) - iVSf{t, x)-2a- VS{t, x) 

Conversely given (0, (j)) as in (B4) satisfying 1)3. 2() it follows that ip satisfies the SE (B2) 
with V as in Q (note R = (l/2)/oc/((^(/>) and S = (l/2)/oc/((/>/0) with exp{R) = {4>(p)^/^).m 

We will discuss this later in more detail and give proofs along with probabilistic content 
(note that the equations (|3.2j) are not imaginary time SE). From this one can conclude 
that nonrelativistic QM is diffusion theory in terms of Schrodinger processes (described 
by {4>,<l)) - more details later). Further it is shown that key postulates in Nelson's sto- 
chastic mechanics or Zambrini's Euclidean QM (cf. |144j ) can both be avoided in con- 
necting the SE to diffusion processes (since they are automatically valid). Look now at 
Theorem 3.1 for one dimension and write T = ht with X = {h/y/m.)x; then the SE 
(B2) becomes (B5) ihipT = — (^^/2?tt.)V'xx — iAtpx + V^P where A = ah/^/m. In ad- 
dition (B6) ihRr + {h^/m'^)RxSx + {h^/2m^)Sxx + ARx = and (B7)F = -ihSr + 
{h^/2m)Rxx + ih^/2m^)Rl - {h^ /2m'^)S\ - ASx- Hence 

PROPOSITION 3.1. Equation (B2), written in the variables (B8) X = {h/y/m)x, T = 
ht, with A = {^/m/h)a and V = V{X,T) ~ V{x,t) is equivalent to (B5). 

Making a change of variables in (|3.2|) now, as in Proposition 3.1, yields 
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COROLLARY 3.1. Equation ()3.2|) . written in the variables of Proposition 3.1, becomes 



(3.4) 



-h(j)T + 7:—4>xx 

zm 



A4>x + c4> = 0; 



-V{X, T) - 2hST Sj^ - 2ASx 



m 



Thus the diffusion processes pick up factors of h and hj \/m. ■ 

REMARK 3.1. We extract here from the Appendix to [nUj for some remarks on com- 
peting points of view regarding diffusion and the the SE. First some work of Fenyes |49| is 
cited where a Lagrangian is taken as 

2" 



(3.5) 



L(t) 



dS 1 



2 \2 



where fJ-tix) = exp{2R(t, x)) denotes the distribution density of a diffusion process and V 
is a potential function. The term (B9) Il{fJ,) = (l/2)[(l/2)(V/u/^)]^ is called a diffusion 
pressure and since (l/2)(V/i/;u) VR the Lagrangian can be written as 



(3.6) 



L 



fidx 



Applying the variational principle 5 L{t)dt = one arrives at 



(3.7) ^ + -[{V{R + S)Y 



1 



2 /i 



1 



{V{R + S))-[-^] + {-^] --^ + V = 



2 fj. 



1 A/i 



4 ^ 



which is called a motion equation of probability densities. From this he shows that the 
function ip = exp{R + iS) satisfies the SE (BIO) idt + (1/2)AV' - V{t,x)'4> = 0. Indeed 
putting (B9) and the formula (Bll) {l/2){Aij,/ ij,) + {l/2)AR+{\/ R)^ into ((SSI) one obtains 



(3.8) ^ + ^{VSf - i(Vi?)2 -Ur + V = 

which goes along with the duality relation (B12) Rt + (1/2) AS + VS-VR+b-VR = where 
(B13) u = (l/2)(a+a) = VRandv = (l/2)(a— a) = VS* as derived in the Nagasawa theory. 
Hence ip = exp{R + iS) satisfies the SE by previous calculations. One can see however that 
the equation 1)3. 6() is not needed since the SE and diffusion equations are equivalent and 
in fact the equations of motion are the diffusion equations. Moreover it is shown in [Unj 
that 1)3. 6(1 is an automatic consequence in diffusion theory with V = —c — 2St — (V5)^ and 
therefore it need not be postulated or derived by other means. This is a simple calculation 
from the theory developed above. ■ 

REMARK 3.2. Nelson's important work in stochastic mechanics jlllj produced the SE 
from diffusion theory but involved a stochastic Newtonian equation which is shown in [30] 
to be automatically true. Thus Nelson worked in a general context which for our purposes 
here can be considered in the context of Brownian motions 



(3.9) B{t) = dt + (1/2)A + 6 • V + a • V; B{t) = -dt + (1/2)A - 6 • V + a • V 
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and used a mean acceleration (B14) a{t,x) = —{l/2)[B{t)B{t)x + B{t)B(t)x]. Assuming 
the duality relations (B12) - (B13) he obtains a formula 

(3.10) a{t,x) = -^[B{t){-b + a)+B{b + a)] = bt + {l/2)V{bf - {b + v) x cur 1(b)- 

-[-vt + (l/2)An + (l/2)(a • V)a + (l/2)(a • V)a - (6 • V)t; -{vV)b-vx curl{b)] 
Then it is shown that the SE can be deduced from the stochastic Newton's equation 

07 -I 

(3.11) a{t, x) = -VV + jr + 7TV(fe^) - {b + v) x curl{b) 

Nagasawa shows that this serves only to reproduce a known formula for V yielding the 
SE; he also shows that (|3.1fl|) also is an automatic consequence of the duality formulation 
of diffusion equations above. This equation (|3.1fl|) is often called stochastic quantization 
since it leads to the SE and it is in fact correct with the V specified there. However the 
SE is more properly considered as following directly from the diffusion equations in duality 
and is not correctly an equation of motion. There is another discussion of Euclidean QM 
developed by Zambrini HSj. This involves (B15) a{t,x) = {l/2)[B{t)B{t)x + B{t)B(t)x] 
(with {aa^y^ = 6'^). It is postulated that this equals (B16) -Vc + bt + {l/2)V{bf - 
b + v) X curl{b) which in fact leads to the same equation for V as above with V = —c — 
2St - (VS')2 ~2b-VS so there is nothing new. Indeed it is shown in j^O] that (B16) holds 
automatically as a simple consequence of time reversal of diffusion processes. ■ 



3.1. SCALE RELATIVITY. There are several excellent and exciting approaches here. 
The method of Nottale [UHl is preeminent (cf. also lEOl IISll Il22]) and 

there is also a nice derivation of a nonlinear SE via fractal considerations in [22] (indicated 
below). The most elaborate and rigorous approach is due to Cresson '33', with elaboration 
and updating in ^ ^ We refer here to [TH [HI 1211 ^ Hi lUi- There are 
various derivations of the SE and we follow |114j here (cf . also |115| I135j ) . The philosophy 
is discussed in [T3| [T H l33 l l3 Hlll4j and we just write down equations here. First a bivelocity 
structure is defined (recall that one is dealing with fractal paths). One defines first 

(3.12) ^^(^)-^— 

d- / y{t) - yjt - At) \ 
-y{t) = hm^,^o, [ ) 

Applied to the position vector x this yields forward and backward mean velocities, namely 
(B17) {d+/dt)x{t) = and {d^/dt)x{t) = 6_. Here these velocities are defined as the 
average at a point q and time t of the respective velocities of the outgoing and incom- 
ing fractal trajectories; in stochastic QM this corresponds to an average on the quantum 
state. The position vector x{t) is thus "assimilated" to a stochastic process which satis- 
fies respectively after (dt > 0) and before {dt < 0) the instant t a relation (B18) dx{t) = 
b+[x{t)]dt + (i^+(t) = b^[x{t)]dt + d^-{t) where C(i) is a Wiener process (cf. [111]). It is 
in the description of ^ that the D = 2 fractal character of trajectories is inserted; indeed 
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that 1^ is a Wiener process means that the d^'s are assumed to be Gaussian with mean 0, 
mutually independent, and such that 

(3.13) < d^+i{t)d^+j{t) >= 2V5ijdt; < d^^i{t)d^^j{t) >= -2VSijdt 

where < > denotes averaging and D is the diffusion coefficient. Nelson's postulate (cf. 
[Ill ) is that T> = h/2m and this has considerable justification (cf. [114j ). Note also that 
(|3.13j) is indeed a consequence of fractal (Hausdorff ) dimension 2 of trajectories follows from 
< > /dt'^ = dt~^, i.e. precisely Feynman's result < v"^ y'^/'^r^ 6t~^/'^ (the discussion 
here in jll4j is unclear however - cf. P). Note also that Brownian motion (used in Nelson's 
postulate) is known to be of fractal (Hausdorff) dimension 2. Note also that any value of 
D may lead to QM and for T> ^ the theory becomes equivalent to the Bohm theory. 
Now expand any function f{x,t) in a Taylor series up to order 2, take averages, and use 
properties of the Wiener process ^ to get 

(3.14) ^ = ^dt + b+-V + VA)f; ^ = (5i + 6_.V-PA)/ 

Let p{x,t) be the probability density of x{t); it is known that for any Markov (hence 
Wiener) process one has (B19) dtp + div[pbj^) = VIS.p (forward equation) and (B20) dtp + 
div{pb-) = —PA/9 (backward equation). These are called Fokker-Planck equations and one 
defines two new average velocities (B21) V = (l/2)[6-|_ + and U = (1/2) [6+ — Con- 
sequently adding and subtracting one obtains (B22) pt + div{pV) = (continuity equation) 
and (B23) div{pU) - VAp = which is equivalent to (B24) div[p{U - Wlog{p))] = 0. 
One can show, using (|3.14^) that the term in square brackets in (B24) is zero leading to 
(B25) U = Wlog{p). Now place oneself in the (C/, V) plane and write (B26) V = V -iU. 
Then write (B27) { dy/dt ) = {l/2){d++d-)/dt and {du/dt) = {l/2){d+ - d-)/dt. Combin- 
ing the equations in one defines (B28) {dv/dt) = dt + V-V and {du/dt) =VA + U-V; 
then define a complex operator (B29) {d' /dt) = {dy/dt) — i{dii/dt) which becomes 

(3.15) i=(i-^VA\+V■V 



dt \dt J 

One now postulates that the passage from classical mechanics to a new nondifferentiable 
process considered here can be implemented by the unique prescription of replacing the 



standard d/dt by d' /dt. Thus consider (B30) S = (^f^^ C{x,V,t)dt^ yielding by least ac- 
tion (B31) {d'/dt){^C/^V^) = dC/dxi. Define then Vi = d£/dVi leading to (B32) V = 
V5 (note this is S and not S). Now for Newtonian mechanics write (B33) L{x,v,t) = 
{l/2)mv^ - U which becomes C{x,V,t) = (l/2)mV^ - il leading to (B34) - Vil = 
m{d' /dt)V. One separates real and imaginary parts of the complex acceleration 7 = {d'V/dt 
to get 

(3.16) d'V = {dv - idu){V - iU) = {dyV - dyU) - i{duV + dyU) 

The force F = —Vil is real so the imaginary part of the complex acceleration vanishes; 
hence 

(3.17) '^V + =^ + U -VV + V -VU + VAV = 
dt dt dt 
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from which dU /dt may be obtained. Differentiating the expression U = Wlog{p) and using 
the continuity equation yields another expression (B35) {dU /dt) = —W{divV) — V{y-U). 
Comparison of these relations yields ^{divV) = Al/ — U A curlV where the curlU term 
vanishes since U is a gradient. However in the Newtonian case V = mV so (B32) implies 
that V is a gradient and hence a generalization of the classical action S can be defined via 
(B36) V = 2WS (note then V{divV) = and curlV = 0). Combining this with the ex- 
pression for U one obtains (B37) S = log{p^^^) + iS. One notes that this is compatible with 
jlllj for example. The way to the SE is now short; set (B38) ip = ^exp{iS) = exp{iS) 
with (B39) V = -liVVilogip) (note U = VVlog{p), V = 2VVS, V = -liWlogip = 
-WVlog^p) + 2PV5 = V -iU); thus for V = mV the relation (B40) V ~ -ihV or Vip = 
—iKVip has a natural interpretation. Putting ip in (B34), which generalizes Newton's law 
to fractal space the equation of motion takes the form (B41) VH = 2iT>m{d' / dt){Vlog{ip)). 
Noting that d' and V do not commute one replaces d' /dt by ()3.15|) to obtain 

(3.18) Vil = 2iVm [dtVlogitp) - iVA{Vlog(ip)) - 2iV{Vlog{ip) ■ V){Vlogi'tp)] 
This expression can be simplified via 

(3.19) VA = AV; (V/ • V)(V/) = (l/2)V(V/)2; ^ = Alogif) + {Vlog{f)f 
This implies 

(3.20) \Aiyiog{iP)) + (VlogiiP) ■ V){Vlog{iP)) = 

Integrating this equation yields (B42) V'^Atp + iDdtip — (pd/2m)'p = up to an arbitrary 
phase factor a{t) which can be set equal to by a suitable choice of phase S. Replacing 
P by h/2m one arrives at the SE (B43) ihipt = —{fi^/2m)A%p + This suggests an 
interpretation of QM as mechanics in a nondifferentiable (fractal) space. 

REMARK 3.3. Some of the relevant equations for dimension one are collected together 
in Section 6. We note that it is the presence of it derivatives that makes possible the intro- 
duction of a complex plane to describe velocities and hence QM; one can think of this as 
the motivation for a complex valued wave function and the nature of the SE. ■ 



We go now to [23] and will sketch some of the material. Here one extends ideas of 
Nottale and Ord in order to derive a nonlinear Schrodinger equation (NLSE). Using the 
hydrodynamic model in |124j one added a hydrostatic pressure term to the Euler-Lagrange 
equations and another possibility is to add instead a kinematic pressure term. The hy- 
drostatic pressure is based on an Euler equation —Vp = pg where p is density and g the 
gravitational acceleration (note this gives p = pgx in 1-D). In |124j one took p = 'ip*ip, h a 
mass-energy parameter, and p = p; then the hydrostatic potential is (for pQ = 1) 

(3.21) b j g{x) ■dr = -b j ^ ■ dr = -blog{p/po) = -hlog{iP*iP) 

Here —blog{'p*%p) has energy units and explains the nonlinear term of |H1 which involved 

r)ih h"^ 

(3.22) in^ = -_^^^ + UiP-b[log{m^ 
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A derivation of this equation from the Nelson stochastic QM was given by Lemos (cf. j79j ) . 
There are however some problems since this equation does not obey the homogeneity con- 
dition saying that the state X\ip > is equivalent to \ip >; however H3.22() is not invariant 
under ip Xtp. Further, plane wave solutions to ()3.22|) do not seem to have a physi- 
cal interpretion due to extraneous dispersion relations. Finally one would like to have a 
SE in terms of alone. Note that another NLSE could be obtained by adding kinetic 
pressure terms (l/2)/?t>^ and taking p = a'4)*il) where v = p/m. Now using the relations 
from HJ theory (B44) (tp/ip*) = exp[2i&(x)/h] and p = 'V&{x) = mv one can write 
(B45) V = —i{h/2m)Vlog{'4>/4'*) so that the energy density becomes (B46) (l/2)p|t;p = 
[afi? /Sm?)ilj'ilj*Vlog{ilj /il)*) ■ Vlog{ilj* /iIj). This leads to a corresponding nonlinear potential 
associated with the kinematical pressure via (B47) (ah? /8m?)Vlog{tp/ip*) • Vlog{'4j* /iIj). 
Hence a candidate NLSE is 



(3.23) ihdt = - — V> + C/^ - bilogim-^ + ^ [^^"9^ ' ^^""9^ 

(apparently this equation has not yet been derived in the literature). Here the Hamiltonian 
is Hermitian and a ^ h are both mass-energy parameters to be determined experimen- 
tally. The new term can also be written in the form (B48) ^log[ip/'il)*) ■ ^log{il}* = 
— \Vlog{iJj/4'*y\^- The goal now is to derive a NLSE directly from fractal space time dynam- 
ics for a particle undergoing Brownian motion. This does not require a quantum potential, 
a hydrodynamic model, or any pressure terms as above. 

REMARK 3.4. One should make some comments about the kinematic pressure 
terms (B49) {l/2)pv'^ ^=> {h^ /2m){a/m)\Vlog{il))\'^ versus hydrostatic pressure terms of 
the form (B50) jCVp/p) <^=^> —hlog[ij)*ip). The hydrostatic term breaks homogeneity 
whereas the kinematic pressure term preserves homogeneity (scaling with a A factor). The 
hydrostatic pressure term is also not compatible with the motion kinematics of a particle 
executing a fractal Brownian motion. The fractal formulation will enable one to relate the 
parameters a, b to h. ■ 

Following Nottale nondifferentiability implies a loss of causality and one is thinking 
of Feynmann paths with < >oc (dx/dt)"^ oc dt"^^^^^^^^^^ with D = 2. Now a frac- 
tal function f{x,e) could have a derivative df/de and renormalization group arguments 
lead to (B51) {df{x,e)/dloge) = a{x) + bf{x,e) (cf. jll4) '). This can be integrated to 
give (B52) f{x,e) = /o(x)[l — C{x){X/e)^^]. Here A~^C(^) is an integration constant and 
/o(x) = —a{x)/b. This says that any fractal function can be approximated by the sum 
of two terms, one independent of the resolution and the other resolution dependent; one 
expects C{x) to be a flucuating function with zero mean. Provided a ^ and b < one 
has two interesting cases (i) e << A with f{x,e) ~ fo{x){X/e)~'' and (ii) e >> A with f 
independent of scale. Here A is the deBroglie wavelength. Now one writes 

(3.24) r{t + dt,dt) -r{t,dt) = b+{r,t)dt + ^+{t,dt) (—^ ; 



r{t,dt) - r{t - dt,dt) - b-{r,t)dt + C-{t,dt) ( 



To/ 
dt\ 
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where (3 = 1/D and 6-|- are average forward and backward velocities. This leads to 
(B53) v±{r,t,dt) = b±{r,t) + ^±{t, dt){dt / tqY~^ . In the quantum case D = 2 one has /? = 
1/2 so dt^~^ is a divergent quantity (so nondifferentiability ensues). Following [7^1 11141 llllj 
one defines 

.oo.A d±r{t) /r(t + At)-r(t) 

(3.25) — — — = limAt-*±o 



(3.27) -VU = -2im[DdtVlog{7p)] - 2DV D 



dt "^^^ \ At 

from which (B54) d±r{t)/dt = b±. Now following Nottale one writes 

rs 26) — = l(^ + —\ -if^-'L. 

^ ' ' dt 2 ydt dt J 2 \ dt dt 

which leads to (B55) {6/dt) = (d/dt) + v - V — iDV^. Here in principle D is a real valued 
diffusion constant to be related to h. (A symbol D for the fractal dimension is no longer 
needed here (?) - e.g. D = 2 with (B56) < dC±idC±j >= ±2D5ijdt.) Now for the 
complex time dependent wave function we take ip = exp[i&/2mD] with p = VS so that 
(B57) V = —2iD^log{ip). The SE is obtained from the Newton equation [F = ma) via 
(B58) - VU = m{5/dt)v = -2imD{5/dt)Vlog{ip). Inserting (B55) gives 

(see |114j for identities involving V). Integrating 1)3. 27() yields (B59) D'^V'^ip + iDdtip — 
{U/2m)'ip = up to an arbitrary phase factor which may be set equal to zero. Now replacing 
D by h/2m one gets the SE (B60) ihdtip + {fi^ /2m)V'^ip = Uij^. Here the Hamiltonian is 
Hermitian, the equation is linear, and the equation is homogeneous of degree 1 under the 
substitution ip — > A^. 

Next one generalizes this by relaxing the assumption that the diffusion coefficient is 
real. Some comments on complex energies are needed - in particular constraints are often 
needed (cf. |129j ). However complex energies are not alien in ordinary QM (cf. !23j 
for references). Now the imaginary part of the linear SE yields the continuity equation 
dtp + V • (pv) = and with a complex potential the imaginary part of the potential will 
act as a source term in the continuity equation. Instead of (B61) < dC±dC± >= ziz2Ddt 
with D and 2mD = h real one sets (B62) < dC±dC± >= ±{D + D*)dt with D and 
2mZ) = h = a+ip complex. The complex time derivative operator becomes (B63) {5/dt) = 
dt + v-V- {i/2){D + Z)*)V^. Writing again (B64) ^ = exp[i&/2mD] = exp{i&/h) one 
obtains (B65) v = —2iDV log[ijj) . The NLSE is then obtained (via the Newton law) as 
(B66) - VU = m{5/dt)v = -2imD{6/dt)Vlog{'ip). Inserting (B63) one gets 
(3.28) 



VU = 2im 



DdtVlog{tp) - 2iD^{Vlog{'^) ■ V){Vlog{ip) - -{D + D*)DV'^{Vlog{il;) 



Now using the identities (i) VV^ = V^V, (ii) 2{Vlog{'4)) ■ V)(yiog{'ip) = ViVlogiip))^ and 
(iii) V'^log{ip) = V'^ip/ip — {Vlog{'il:))'^ leads to a NLSE with nonlinear (kinematic pressure) 
potential, namely 

(3.29) ihdtil^ = - — + f/V' - i—^Uviogm^i; 

2m n 2m n 



18 



ROBERT CARROLL UNIVERSITY OF ILLINOIS, URBANA, IL 61801 



Note the crucial minus sign in front of the kinematic pressure term and also that h = 
a + i(3 = 2mD is complex. When /3 = one recovers the linear SE. The nonlinear potential 
is complex and one defines (B67) W = —{h? /2m){l3/h){Vlog{7p))'^ with U the ordinary 
potential; then the NLSE is (B68) ihdtij = [-{h^ /2m){a/h)V'^ + U + iW]^). This is the 
fundamental result of it has the form of an ordinary SE with complex potential U + iW 
and complex h. The Hamiltonian is no longer Hermitian and the potential itself depends 
on ip. Nevertheless one can have meaningful physical solutions with real valued energies 
and momenta; the homogeneity breaking hydrostatic pressure term —b{log{ip*'il))ip is not 
present (it would be meaningless) and the NLSE is invariant under ip — > Xtp. 

REMARK 3.5. One could ask why not simply propose as a valid NLSE an equation 

(3.30) ihdti> = -I^VV + Ui, + ^-\Vlog{i;)\^^ 

2m 2m m 

Here one has a real Hamiltonian satisfying the homogeneity condition and the equation 
admits soliton solutions of the form (B69) = CA{x — vt)exp[i{kx — Lot)] where A{x — vt) 
is to be determined by solving the NLSE. The problem here is that the equation suf- 
fers from an extraneous dispersion relation. Thus putting in the plane wave solution 
"0 ~ exp[—i{Et — px)] one gets an extraneous EM relation (after setting U = 0), namely 
(B70) E = (p^/2m)[l + (a/m)] instead of the usual E = p^/2m and hence Eqm / EpT 
where FT means field theory. ■ 

REMARK 3.6. It has been known since e.g. |129j that the expression for the energy 
functional in nonlinear QM does not coincide with the QM energy functional, nor is it 
unique. To see this write down the NLSE of [Hj in the form (B71) ihdtip = dH{ip, 'ilj*)/d'ip* 
where the real Hamiltonian density is 

(3.31) Hiip, V*) = ^A*VV + Uip*^; - b^*log('ip*tp)ip + bip^ip 

2m 

Then using Epx = f Hd^r we see it is different from < H >qm and in fact EpT — Eqm = 
J hil)*iljd?r = b. This problem does not occur in the fractal based NLSE since it is written 
entirely in terms of V'- B 

REMARK 3.7. In the fractal based NLSE there is no discrepancy between the QM 
energy functional and the FT energy functional. Both are given by 

(3.32) Nf,.^ji, = -|!-^^*vv + uri^ - i^^riyiogii^fi' 

The NLSE is unambiguously given by (B71) and H{'iIj,iJj*) is homogeneous of degree 1 in 
A. Such equations admit plane wave solutions with dispersion relation E = p^ /2m; indeed, 
inserting the plane wave solution into the fractal based NLSE one gets (after setting [7 = 0) 

(3 33) E- ^ I _ p^ a + il3 _ p^ 

2m h 2m h 2m 2m h 2m 
since h = a + if3. The remarkable feature of the fractal approach versus all other NLSE 
considered sofar is that the QM energy functional is precisely the FT one. The complex 
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diffusion constant represents a truly new physical phenomenon insofar as a small imaginary 
correction to the Planck constant is the hallmark of nonlinearity in QM (see j23i for more 
on this). ■ 

4. REMARKS ON A FRACTAL SPACETIME 

There have been a number of articles and books involving fractal methods in spacetime or 
fractal spacetime itself with impetus coming from quantum physics and relativity. We refer 
here especially to P III L3, ,24, .58, ,94, 95, 96, 97, 98, 99, 100, lOHHnilLOS, lOSjIMl for 
background to this paper. Many related papers are omitted here and we refer in particular 
to the journal Chaos, Solitons, and Fractals CSF) for further information. For information 
on fractals and stochastic processes we refer for example to [U H?! |2H1 HH ISl IHH IZH UHl 
IHa ESI HUB Uni ins inn Unni CSni IHSl • We discuss here a few background ideas and 
constructions in order to indicate the ingredients for El Naschie's Cantorian spacetime 
whose exact nature is elusive. Suitable references are given but there are many more papers 
in the journal CSF by El Naschie (and others) based on these fundamental ideas and these 
are either important in a revolutionary sense or a fascinating refined form of science fiction. 
In what appears at times to be pure numerology one manages to (rather hastily) produce 
amazingly close numerical approximations to virtually all the fundamental constants of 
physics (including string theory). The key concepts revolve around the famous golden ratio 
("v/S — l)/2 and a strange Cantorian space which we try to describe below. It is very 
tempting to want all of these (heuristic) results to be true and the approach seems close 
enough and universal enough to compel one to think something very important must be 
involved. Moreover such scope and accuracy cannot be ignored so we try to examine some 
of the constructions in a didactic manner in order to possibly generate some understanding. 

4.1. COMMENTS ON CANTOR SETS. 

EXAMPLE 4.1. In the paper |^j one discusses random recursive constructions leading 
to Cantor sets, etc. Associated with each such construction is a universal number a such 
that almost surely the random object has Hausdorff dimension a (we assume that ideas of 
Hausdorff and Minkowski-Bouligand (MB) or upper box dimension are known - cf. 1141 
1461 180j ) . One construction of a Cantor set goes as follows. Choose x from [0, 1] according to 
the uniform distribution and then choose y from [x, 1] according to the uniform distribution 
on [x, 1]. Set Jo = [0, x] and Ji = [y, 1] and recall the standard 1/3 construction for Cantor 
sets. Continue this procedure by rescaling to each of the intervals already obtained. With 
probability one one then obtains a Cantor set with Hausdorff dimension (CI) a = (p = 
("v/S — l)/2 ~ .618. Note that this is just a particular random Cantor set; there are others 
with different Hausdorff dimensions (there seems to be some - possibly harmless - confusion 
on this point in the El Naschie papers). However the golden ratio is a very interesting 
number whose importance rivals that of vr or e. In particular (cf. |3) <^ is the hardest number 
to approximate by rational numbers and could be called the most irrational number. This 
is because its continued fraction represention involves all I's. ■ 

EXAMPLE 4.2. From the Hausdorff (H) dimension of a traditional triadic Cantor 

set is di^'* = log{2)/log{3). To determine the equivalent to a triadic Cantor set in 2 di- 
mensions one looks for a set which is triadic Cantorian in all directions. The analogue 
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of an area A = 1 x 1 is a quasi-area Ac = dc x dc and to normalize Ac one uses 
P2 = {A/Ac)2 = (for n-dimensions (C2) /)„ = Then the n*'^ Can- 

tor hke H dimension dc"^ will have the form (C3) di"^ = Pnd^^ = l/(dc'^^)"~^. Note also 
that the H dimension of a Sierpinski gasket is (C4) di"^^V^c"^ = ^/d^c^ = log{2>) /log{2) 
and in any event the straight-forward interpretation of dc = log{3)/log{2) is a scal- 
ing of d^c^ = log{2)/log(3) proportional to the ratio of areas {A/Ac)2- One notes that 
(C5) di^^ = l/(di°V = ilog{3)/log{2)f ~ 3.997 ~ 4 so the 4-dimensional Cantor set is 
essentially "space fihing". 

Another derivation goes as follows. Define probability quotients $7 = dim{subset) /dim[set). 
For a triadic Cantor set in 1-D (C6) = d^c'' /d^P = d^°^ {d^P = 1). To hft the Cantor set 
to n-dimensions look at the multiplicative probability law (C7) O*^"-* = (fi'-"'^-*)" = {d^c^)"^. 
However since = df^ /d^c^ we get (C8) df^ /d'^c^ = (di°^)" ^ di") = l/(di°^)"-i. Since 
j^(n-i) jg probability of finding a Cantor point (Cantorian) one can think of the H di- 
mension di"^ = l/r^^""^) as a measure of ignorance. One notes here also that for di'^^ = (f) 

(the Cantor set Sc'^^ of Example 2.1) one has d^f^ = l/(/)'^ = 4-|-0^~ 4.236 which is surely 
space filling. ■ 



Based on these ideas one proves in |951 1961 1^5] a number of theorems and we sketch some of 

this here. One picks a "backbone" Cantor set with H dimension d^c^ (the choice of (/> = d^c"^ 
will turn out to be optimal for many arguments). Then one imagines a Cantorian spacetime 
built up of an infinite number of spaces of dimension di"^ (— oo < n < oo). The exact 
form of embedding etc. here is not specified so one imagines e.g. 

unions and intersections) in some amorphous sense. There are some connections of this to 
vonNeumann's continuous geometries indicated in In this connection we remark that 

only G;^"""-* is the completely empty set (G;^""'^^ is not empty). First we note that (f)'^+(j)—l = 
leading to (C9) 1 + (P= 1/4), 0^ = (2 + ,/.)/,^, (1 + </.)/(!-(/)) = l/(/)(l - ,^) = 4 + </>3 = 1/4)^ 
(a very interesting number indeed). Then one asserts that 



THEOREM 4.1. Let {Vl'-^^Y be a geometrical measure in n-dimensional space of a multi- 
plicative point set process and fi^^-* be the Hausdorff dimension of the backbone (generating) 
set di°\ Then < d >= l/di°^(l - df^) (called curiously an average Hausdorff dimension) 
will be exactly equal to the average space dimension - < n >= {1 + di'^^)(l — d^c"*) and 
equivalent to a 4-dimensional Cantor set with H-dimension dc^^ = l/(dc'^^)^ if and only if 

df^ = 4>. 

To see this take = (0(i))" again and consider the total probability of the additive 
set described by the namely (CIO) Zq = = 1/(1 -S^^^^)- It is conceptually 

easier here to regard this as a sum of weighted dimensions (since d^'^ = l/{d^c^)'^~^) and 
consider Wn = n{d^c^)^. Then the expectation of n becomes (note d^^ ~ l/{d^^)^^^ ~ 
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Another average here is defined via (blackbody gamma distribution) 
^ ' f^n{mYdn log{Q(^)) 

which corresponds to ~ < n > after expanding the logarithm and omitting higher order 
terms. However ~< n > seems to be the more vahd calculation here. Similarly one defines 
(somewhat ambiguously) an expected value for dc"^ via 

(4.3) <d>= ^' ""^^ = 

This is contrived of course (and cannot represent E{d^^) since one is computing reciprocals 
^(n/di"^) but we could think of computing an expected ignorance and identifying this with 
the reciprocal of dimension. Thus the label < d > does not seem to represent an expected 
dimension but if we accept it as a symbol then for dc"^ = </> one has from (C9) 

(4.4) ^<„>=i±|=<.>=^=4«=4 + #» = Jj~4.236 

REMARK 4.1. We note that the normalized probability (Cll) N = Q'-^^/Zq = 
0(i)(l - J^(i)) = 1/ < d > for any Further if < d >= 4 = l/dlP\l - d^P^) one 

has d^'^ = 1/2 while- < n >= 3 < 4 =< d >. One sees also that d^^ = 1/2 is the minimum 
{where d < d> /d{d^c^) =0). U 

REMARK 4.2. The results of Theorem 4.1 should really be phrased in terms of (2;°° 
(cf. |inij ). thus {H ~ Hausdorff dimension and T ~ topological dimension) 

(4.5) ei^^^gW = 4") = ^ ; 

1 1 + d(°) 

In any event is formally infinite dimensional but effectively it is 4 it dimensional with an 
infinite number of internal dimensions. We emphasize that appears to be constructed 
from a fixed backbone Cantor set with H dimension 1/2 < d^c^ < 1; thus each such d^^ 
generates an space. Note that in jlUl] is looked upon as a transfinite discretum (?) 
underpinning the continuum. ■ 

REMARK 4.3. An interesting argument from jlUUj goes as follows. Thinking of d^c^ as a 
geometrical probability one could say that the spatial (3-dimensional) probability of finding 
a Cantorian "point" in must be given by the intersection probability (C12) P = (d^c^)^ 
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where 3 ~ 3 topological spatial dimension. P could then be regarded as a Hurst exponent 
(cf. [H I114| ri43j ) and the Hausdorff dimension of the fractal path of a Cantorian would be 
(C13) dpath = '^/H = l/P = 1/(4°^)^. Given = cj) this means dpath = 4 + ^^ ^ 4+ so a 
Cantorian in 3-D would sweep out a 4-D world sheet; i.e. the time dimension is created by 
the Cantorian space (! - ?). Conjecturing further (wildly) one could say that perhaps 
space (and gravity) is created by the fractality of time. This is a typical form of conjecture 
to be found in the El Naschie papers - extremely thought provoking but ultimately heuristic. 
Regarding the Hurst exponent one recalls that for Feynmann trajectories in 1 + 1 dimensions 
(C14) dpath = ^/H = = d^\ Thus we are concerned with relating (C13) and (C14) 

(among other matters). Note that path dimension is often thought of as a fractal dimension 
(M-B or box dimension), which is not necessarily the same as the Hausdorff dimension. 
However in ^ one shows that quantum mechanical free motion produces fractal paths of 
Hausdorff dimension 2 (cf. also [ZOl)- ^ 

REMARK 4.4. Following |25j let Sc'^^ correspond to the set with dimension d^^ = <j). 

Then the complementary dimension is d^i^ = 1 — cf) = cjp' . The path dimension is gien as in 

(C14) by (C15) dpath = d^c^ = 1/0=1 + and dpath = dP = 1/(1-0) = V^' = (1 + 0)'- 
Following El Naschie for an equivalence between unions and intersections in a given space 
one requires (in the present situation) that 

(4.6) d.rU = d(^^ + S^^='- + ^ = ^^ = ^='-.^ = dp.S^^=i + c^' 

where (C16) dcrit = 4 + 0^ = d^^ ~ 4.236. Thus the critical dimension coincides with the 
Hausdorff dimension of S'c^^ which is embedded densely into a smooth space of topological 
dimension 4. On the other hand the backbone set of dimension d^'^ = (j) is embedded 
densely into a set of topological dimension zero (a point). Thus one thinks in general of 

(n) 

dc as the H dimension of a Cantor set of dimension embedded into a smooth space of 
integer topological dimension n. ■ 

REMARK 4.5. In [2S1 it is also shown that realization of the spaces ^^""^ comprising 
can be expressed via the fractal sprays of Lapidus-van Frankenhuysen (cf. lEHI)- Thus 
we refer to |8f)j for graphics and details and simply sketch some ideas here (with apologies 
to M. Lapidus). A fractal string is a bounded open subset of R which is a disjoint union 
of an infinite number of open intervals £ = ii,£2,---- The geometric zeta function of 
£ is (C17) C£(s) = YIT '^J'^- One assumes a suitable meromorphic extension of C£ and 
the complex dimensions of £ are defined as the poles of this meromorphic extension. The 
spectrum of £ is the sequence of frequencies f = k ■ {k = 1,2, ■ ■ ■ ) and the spectral zeta 
function of £ is defined as (C18) Cu{s) = J2f where in fact G(s) = C£(s)C(s) (with ({s) 
the classical Riemann zeta function). Fractal sprays are higher dimensional generalizations 
of fractal strings. As an example consider the spray obtained by scaling an open square 
B of size 1 by the lengths of the standard triadic Cantor string CS. Thus f2 consists of one 
open square of size 1/3, 2 open squares of size 1/9, 4 open squares of size 1/27, etc. (see jHO] 
for pictures and explanations). Then the spectral zeta function for the Dirichlet Laplacian 
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on the square is (C19) Cb{s) = X^J^J n2=i("'i ^ n'^Y^'^ and the spectral zeta function of the 
spray is (C20) Cu{s) = Ccs{s) • Cb{s)- Now 6;°° is composed of an infinite hierarchy of sets 
<£^^^ with dimension (1 + (/')-'~^ = l/(f)^~^ {j = 0, ±1, ±2, • • • ) and these sets correspond to a 
special case of boundaries for fractal sprays whose scaling ratios are suitable binary 
powers of 2~'f''' . Indeed for n = 2 the spectral zeta function of the fractal golden spray 
indicated above is (C21) Ci^(s) = (1/(1 — 2 • 2*'^)^s(s). The poles of Cb{s) do not coincide 
with the zeros of the denominator 1 — 2 • 2~''^'^ so the (complex) dimensions of the spray 
correspond to those of the boundary dO, of 0,. One finds that the real part 5Rs of the complex 
dimensions coincides with dim(E^'^^ = 1 + (j) = l/cfP' and one identifies then 517 with 
The procedure generalizes to higher dimensions (with some stipulations) and for dimension 
n there results SRs = l/<j)'^~^ = dim<£.^"'\ This produces a physical model of the Cantorian 
fractal space from the boundaries of fractal sprays (see [2^] for further details and [HOI for 
precision). Other (putative) geometric realizations of <t°° are indicated in |in4j in terms of 
wild topologies, etc. ■ 

5. HYDRODYNAMICS AND THE FRACTAL SCHRODINGER EQUATION 

We sketch first some material from [Hj (see also |141 1114| I115( I116j and Sections 2-4 for 
background). Thus let ip be the wave function of a test particle of mass tuq in a force 
field U{r,t) determined via (Dl) ihdtip = Uif) — (/i^/2m)V^'0 where = A. One writes 
(D2) Tp{r,t) = R{r,t)exp{iS{r,t)) with (D3) v = {h/2m)VS and p = R - R (one assumes 
p ^ for physical meaning). Thus the field equations of QM in the hydrodynamic picture 
are 

(5.1) dt{mopv) = dt{mopv) + \/{mopv) = -pV{U + Q); dtp + V ■ {pv) = 

where (D4) Q = — (/i^/2mo)(A-^/p/y^) is the quantum potential (or interior potential). 
Now because of the nondifferentiability of spacetime an infinity of geodesies will exist be- 
tween any couple of points A and B. The ensemble will define the probability amplitude 
(this is a nice assumption but what is a geodesic here). At each intermediate point C one 
can consider the family of incoming (backward) and outgoing (forward) geodesies and define 
average velocities 6+(C) and 6_(C) on these families. These will be different in general and 
following Nottale this doubling of the velocity vector is at the origin of the complex nature 
of QM. Even though Nottale reformulates Nelson's stochastic QM the former's interpre- 
tation is profoundly different. While Nelson (cf. |lllj ^ assumes an underlying Brownian 
motion of unknown origin which acts on particles in a still Minkowskian spacetime, and 
then introduces nondifferentiability as a byproduct of this hypothesis, Nottale assumes as 
a fundamental and universal principle that spacetime itself is no longer Minkowskian nor 
differentiable. While with Nelson's Browian motion hypothesis, nondifferentiability is but 
an approximation which expected to break down at the scale of the underlying collisions 
(?), where a new physics should be introduced, Nottale's hypothesis of nondifferentiability 
is essential and should hold down to the smallest possible length scales. (This sentence is 
interesting but needs elaboration). Following Nelson one defines now the mean forward and 
backward derivatives 

/.o^ d± / y{t + At) - y{t) 

(5.2) —yW = hmAt^o± 
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This gives forward and backward mean velocities (D5) {dj^ / dt)x{t) = and {d^ / dt)x{t) = 
b- for a position vector x. Now in Nelson's stochastic mechanics one writes two systems 
of equations for the forward and backward processes and combines them in the end in a 
complex equation, Nottale works from the beginning with a complex derivative operator 

2^ S _ {d+ + - i{d+ - d_) 



dt 2dt 

leading to (D6) V = {5/dt)x{t) = v - iu = (l/2)(6+ + 6_) - ii/2){b+ - One defines 
also (D7) (d^/dt) = {l/2){d+ + d^)/dt and {dj/dt) = {l/2){d+ - dJ) /dt so that d^x/dt = v 
and dux/dt = u. Here v generalizes the classical velocity while u is a new quantity aris- 
ing from nondifferentiability. This leads to a stochastic process satisfying (respectively for 
the forward {dt > 0) and backward {dt < 0) processes) (D8) dx{t) = b^[x{t)] = 
b-[x{t)] + d£,-{t). The terms can be seen as fractal functions and they amount to a 

Wiener process when D = 2 (presumably the fractal dimension). Then the d(^{t) are Gauss- 
ian with mean zero, mutually independent, and satisfy (D9) < dS^±id£,±j >= ziz2D5ijdt 
where D is a diffusion coefficient. D can be found via D = h/2mQ given tq = h/{moc^) 
(deBroglie time scale in the rest frame - cf ^1] for more on this). Now (D9) allows one to 
give a general expression for the complex time derivative, namely 

df I d'^f 

(5.4) df = ^ + Vf-dx + -—J—dx^dxj 

ot 2 oxiOXj 

Next compute the forward and backward derivatives of /. Then < dxidxj >— >< dS,±id(^-i-j > 
so the last term in 1)5.4(1 amounts to a Laplacian via (D9) and one obtains (DIO) {d±f /dt) = 
[3^ -|- 6-1- • V lb DA]/. This is an important result. Thus assume the fractal dimension is not 
2 in which case there is no longer a cancellation of the scale dependent terms in ()5.4() and 
instead of DAf one would obtain an explicitly scale dependent behavior Ddt^"^/^^"^ Af . 
In other words the value D = 2 implies that the scale symmetry becomes hidden in the 
operator formalism. Using (DIO) one obtains the complex time derivative operator in the 
form (Dll) {5/dt) = dt + V -V - iDA (cf. (D6) for V). Nottale's prescription is then to 
replace d/dt by 5/dt. In this spirit one can write now (D12) tp = exp{i{& /2mQD)) so that 
(D13) V = —2iD'V{log{'ip)) and then the generalized Newton equation (D14) — VfJ = 
niQ{6/dt)V reduces to the SE. 

Now assume the velocity field from the hydrodynamic model agrees with the real part v 
of the complex velocity V and equate the wave functions from the two models (D12) and 
(D2); one obtains for & = s + ia (D15) s = 2mQDS, D = (/i/2mo), and a = —mQDlog{p). 
Using the definition (D16) V = (l/mo)VS = (l/mo)Vs -|- {i/mo)'Va = v — iu (which 
results via (D6) by putting (D12) into (D13)) we get (D17) v = (l/mo)Vs = 2DVS 
and u = — (l/mo)Vcr = D'Vlog{p). Note that the imaginary part of the complex velocity 
given in (D17) coincides with Nottale. Dividing the time dependent SE (Dl) by 2mo and 
taking the gradient gives (D18) VU/niQ = 2DV[idtlog{ip) + D{Aip/'tp)] where h/2mo has 
been replaced by D. Then consider the identities 

(5.5) AV = VA; (V/ • V)(V/) = {l/2)V{Vff; ^ = Alog{f) + {Vlog{f)f 
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Now the second term in the right of (D18) becomes (D19) V(A'0/^) = A{'Vlog{ip)) + 
2{Vlog{ip) ■ V){Vlog{i>)) so (D18) can be written as (D20) VU = 2iDmo[dtVlog{ij) - 
iD/\{Vlog{il)) - 2iD{Vlog{'4)) ■ V){Vlog{i;))]. One can show that (D20) is nothing but the 
generahzed Newton equation (D14). Now if we replace the complex velocity (D13), taking 
into account (D6) and (D17) we get 

(5.6) -VU = mo{dtiv - iDVlog{p) + [i{v - iDVlog{p) • V](t; - iDVlog{p))- 

-iDA{v-iDVlog{p))} 
Equation ()5.6|) is a complex differential equation and reduces to (using 

(5.7) mo[dtv + (v ■ V)v] = -V (^U - 2moD^^^ ; y | ^ [S^p + V • (pv)] | 

The last equation in (|5.7|) reduces to the continuity equation up to a phase factor a{t) which 
can be set equal to zero (note again that p 7^ is posited). Thus (|5.7|) is nothing but the 
fundamental equations (|5.1j) of the hydrodynamic model. Further combining the imaginary 
part of the complex velocity in (D17) with the quantum potential (D4) and using (|5.5|) one 
gets (D21) Q = — mo-DV-M — (l/2)moM^. Since u arises from nondifferentiability according 
to our nondifferentiable space model of QM it follows that the quantum potential comes 
from the nondifferentiability of the quantum spacetime (very nice but where is from 
the title of [3j - also the x derivatives should be clarified). 

Putting [/ = in the first equation of (|5.7j) . multiplying by p, and taking the second 
equation into account yields 

d d 

(5.8) dtirriQpvk) + j—{mQpuiVk) = -p — 



VP 



(here Vj. ~ Vk seems indicated). Now set (D22) = mopviVk — (Tik along with cTjfc = 
mopD'^{d/dxi){d/dxk){log{p)). Then 1)5. 8|) takes the simple form (D23) dtimopuk) = 
—dUik/dxi. The analogy with classical fluid mechanics works well if one introduces the 
kinematic (D24) p = D/2 and dynamic = {l/2)mQDp viscosities. Then Iljfc defines the 
momentum flux density tensor and Uik the internal stress tensor (D25) aik = r][{dui/ dxk) + 
(duk/dxi)]. From (D22) one can see that the internal stress tensor is build up using the 
quantum potential while the equations (|5.H) or (|5.7|) are nothing but systems of Navier- 
Stokes type for the motion where the quantum potential plays the role of an internal stress 
tensor. In other words the nondifferentiability of the quantum spacetime manifests itself like 
an internal stress tensor. For clarity in understanding (D23) we put this in one dimensional 
form so ()5.8)) becomes 

(5.9) dtimopv) + d^imopv^) = -pd (2moD^^d^^) = pdQ 



VP 

and n = rriQpv^ — a with a = mQpD'^d^log(p). This agrees in the standard formulas (cf. 
dH). Now note = (l/2)p-V2p' and = (l/2)[-(l/2)p-3/2(^/)2 ^ ^-1/2^//] ^^^j^ 

d-'logip) = d{p'/p) = {p"/p) - {p'/pf while 



(5.10) -pd 



2moD^^ (92^) 
y/P 



-2moD^pd 



' ^ -Ip-np'f + p-'^'p"' 



2^p\ 2 
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-2moD^pd 



EL 

2p 



p 



—moD^pd 



pL 

p 



'7) 



One wants to show then that (D23) holds or equivalently —da = (|5.1U|) . Here 



(5.11) -da = -d[mQpD^dHog{p)] = -m^D 



— - - 




Pdi'--^-^] 
\P P ) 



so we want (|5.11|) = (|5.1())1 and this is easily verified. 



6. RECAPITULATION 

We write down now some of the main formulas here (with some unification of notation) 
in order to help provide perspective. The goal is not entirely clear but many questions 
will arise as we go along and at the end. Hopefully we will be able to answer some of the 
questions. 

(1) We write from Section 2 (El) ^ = Rexp{iS/fi) with 



(6.1) 



(6.2) 



St 

For P = R^ and Q 



2m 2m R m 







-{tf/2m){R"/R) this yields 



St 



[S 



r\2 



+ Q + V = 0; Pt + -{PSy 
2m m 



Writing p = mP and p = mx leads to 



(6.3) 



VP 







Along the way one arrived at 1)2. 8(1 and "completed" this with a pressure term 
VF = p~^Vp or F' = (l/i?^)p' to arrive at (E2) mvt + mvv' = -d{V + Q) - F' 
corresponding to a SE (E3) i/i^j = —{h'^/2m)ip" + Vip + Fip. One wants then 
F = F{ip). 

(2) Consider a quantum state corresponding to a "subquantum" statistical ensemble 
governed by classical kinetics in a phase space. One arrives at ■0 = p^^'^exp{i&/h) 
with (E4) ihipt = -{h^/2m)ij^^ + Vtp where e = NS, N = /|VpcZ"x, h = 
Nrj, m = Np, V = NV, and log{ip) = {l/2)log{p) + {i/ri)S. The fields p, S or 
S determine a quantum fluid with (cf. (|2.15j) ) 



(6.4) 



di 1 d'^S 1 di dS 

— H 1 ^ — = 0; 

dt p dx^ pdx dx 



dS_ _ rf_ei^ _!f_(^ 
dt 4/i dx"^ 8p \dx 

which for 'ip = p^^'^exp{i&/h) leads to 



+ 



1 fdS 



2p \dx 



+ V = 



(6.5) 



d^^ 

2m dx"^ 



+ V* 
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(3) The Fisher information connection la Remarks 2.4-2.5 involves a classical ensemble 
with particle mass m moving under a potential V 

(6.6) St + 7^(5')^ + V = 0; Pt + -d{PS')' = 

Im m 

where S is a momentum potential; note that no quantum potential is present but 
this will be added on in the form of a term (l/2m) J dt{AN)'^ in the Lagrangian 
which measures the strength of fluctuations. This can then be specified in terms 
of the probability density P as indicated in Remark 2.4 leading to a SE 1)2. 24() . A 
"neater" approach is given in Remark 2.5 leading in 1-D to 

(6.7) St + ^{Sr + V + - f ^ - ^) = 

Note that Q = -{h^ /2m){R" / R) becomes for R = P^/^ (^5) Q = -{2h^ /2m)[{2P" / P)- 
(P'/P)2] (cf. 1)2.4(1 ). Thus the addition of the Fisher information serves to quantize 
the classical system. 

(4) One defines an information entropy (IE) in Remark 2.6 via (E6) & = — j plog{p)d?x {p - 
I'i/'l^) leading to 

(6.8) 9^ = y + log{p))d{vp) ~ y ^ 

modulo constants involving D ~ h/2m. & is typically not conserved and dtp = 
—V • (vp) {u = D'Vlog{p) with v = —u corresponds to standard Brownian motion 
with d&/dt > 0. Then high IE production corresponds to rapid flattening of the 
probability density. Note here also that 5 ~ ~(2/^^) / pQdx = J dx[{p')'^/p] is a 
functional form of Fisher information. Entropy balance is discussed in Remark 2.8 
and the manner in which Q appears in the hydrodynamical formalism is exhibited 
in 

(5) The Nagasawa theory (based in part on Nelson's work) is very revealing and fasci- 
nating (see [nnil^)- The essense of Theorem 3.1 is that = exp{R + iS) satisfies 
the SE (E7) ii;t + (1/2)V'" + iaip' - 't^V = if and only if 

(6.9) V = -St + ^R" + ^{R'f - ^{S'f -aS; = Rt + ^S" + S'R' + aR' 

Changing variables in (E8) {X = i^j ^pm)x and T = h£) one arrives at (E9) ih^T = 
— {h? /2m)ipxx — iAipx + where A = ahj^fm and 

(6.10) ifiRT + {h^/rr?)RxSx + {h^ /2m^)Sxx + ARx = 0; 

V = -ihSx + {h^/2m)Rxx + {h^ /2m^)R\ - {h^ /2m'^)S\ - ASx 
The diffusion equations then take the form 

(6.11) h<t)T + ^(kxx + A(l)x + c(j) = 0; -Hc^t + i^(pxx - AcPx + ccp = 0; 

2m 2m 

/?2 

c = -V(X, T) - 2hST Si- 2ASx 

m 
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It is now possible to introduce a role for the quantum potential in this theory. Thus 
from if) = exp{R + iS) (with h = m = 1 say) we have ip = p^^'^exp{iS) with p^/^ = 
exp{R) orR = {l/2)loglp). Hence {l/2){p'/p) = R' and R" = {l/2)[{p" / p)-{p' / p)'^\ 
while the quantum potential is Q = (1/2) ((92^1/2/^1/2) ^ -{l/^)[{2p" / p) - {p' / pf] 
(cf. (|2.4|) ). Equation becomes then 

1 [2p" {pr\ I,, .2 „c - c 1 



(6.12) V = -St + -[^-^] - -{S'Y -aS = St + -{S'Y + V + Q + aS = 0; 

8 \ p p-' J 2 2 

pt + pS" + S'p' + ap' = = pt + {pS'y + ap' = 
Thus -2St - {S'f = 2V + 2Q + 2AS and one has 

PROPOSITION 6.1. The creation-annihilation term c in the diffusion equations 
(cf. Theorem 3.1) becomes 

(6.13) c = -V -2St- {S'f - 2aS' = V + 2Q + 2a{S - S') 
where Q is the quantum potential. 

(6) Regarding scale relativity one writes (cf. 1)3. 12() 

(6.14) -y{t) = hm^.^o^ 

and we collect equations in {p = j-fAp) 

(6.15) dx = h+dt + di+ = h^dt + d^^;< df^ >= 2Vdt = - < d^^ > 

(6.16) ^ = (9, + b+d + Vd^)f; ^ = {dt + b.d - Vd^)f 

(6.17) V = ^(6+ + 6_); U = ^(5+ - 6_); pt + d{pV) = 0; U = Vd{log{p)); 

V = V-iU; dv = \{d+ + d_); du = \{d+ - d.) 

(6.18) '^ = dt + vd] ^ = vd^ + ud; ^ = {dt-ivd^) + vd 

dt dt dt 

(6.19) V = 2VdS; S = log{p^^^) + iS; ^ = ^e'^ = e'^; V = -2iVdlog{^) 
For Lagrangian C = (l/2)mV2 - mil one gets a SE 

(6.20) ihilJt = + iiV' 

2m 

coming from Newton's law (ElO) — 3il = —2i'Dm{d' / dt)dlog{ilj) = m{d' /dt)V. 
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(7) The development in Section 3 based on 23 involves thinking of nonlinear QM 
as a fractal Brownian motion with complex diffusion coefficient. We note (ElO) 
corresponds to (B58) and (B55) arises in ()6.18|1 . These give rise to 

(6.21) -VU = -2im[DdtVlog{ip)] - 2DV D- ^ 



Thus putting in a complex diffusion coefficient leads to the NLSE 

h? a It? 3 

(6.22) ihdttl; = -^]-vV + U^- ^^|(V^05(V'))V 

with h = a + il3 = 2mD complex. 
(8) In j3j one writes again ip = Rexp{iS/h) with field equations in the hydro dynamical 
picture 

(6.23) dt{mopv) = dt{mopv) + V{mopv) = -pV{u + Q); dtp + V ■ (pv) = 

where Q = — (?i^/2mo)(Ay^/y^). One works with the Nottale approach as above 
with ~ dv and du ~ du (cf. (|6.18j) ). One assumes that the velocity field 
from the hydro dynamical model agrees with the real part v of the complex velocity 
V = V — iu so (cf. (|6.17j) ) V = (l/mo)Vs ~ 2Dds and u = —{l/mo)Va ~ Ddlog{p) 
where D = h/2mQ. In this context the quantum potential Q = — (/i^/2mo)Ay^/y^ 
becomes (Ell) Q = —mQDVu—{\/2)mQV? ~ —{h/2)du—{\/2)mQV?. Consequently 
Q arises from the fractal derivative and the nondifferentiability of spacetime. Further 
one can relate u (and hence Q) to an internal stress tensor (D25) whereas the v 
equations correspond to systems of Navier-Stokes type. Note here that involves 
a term relating the stress tensor 11 and Q directly. 

7. CONCLUSIONS 

One feature either exhibited or suggested in the examples displayed involves the role 
of a quantum potential in either quantization or "classicalization" of certain systems of 
equations of hydrodynamic type. Now with numbers referring to Section 6 we have: 

(1) One arrived at an equation of hydrodynamic type directly from the SE upon addition 
of a pressure term which served to augment the original potential V (however this 
could have simply been included in V). On the other hand Q does not appear in the 
SE but is generated by the decomposition -0 = Rexp{iS/h) 

(2) In a general statistical mechanical approach, with the dynamics determined by clas- 
sical kinetics in a phase space, the quantum potential has an interpretation in terms 
of an internal stress tensor for a quantum fluid. The equations are again described 
in terms of a probability density p and a phase factor S. 

(3) In ^3-#4 one takes a classical statistical ensemble with S a momentum potential 
and expresses momentum fluctuations in terms of Fisher information; this leads to 
a SE with quantization term Q expressed as Fisher information. In Remarks 2.6-2.8 
we show how Fisher information, entropy, and the quantum potential are mutually 
entangled (cf. also [SHI)- In (|2.38|) - ()2.39|) (based on [S5]) we see how the Euler 
equation {dt + v ■ V)v = (F/m) — (AP/p) (where P is a pressure term) is related to 
the quantized form {dt + v ■ V)v = {F/m) — VQ arising from a SE. 
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(4) The Nagasawa-Nelson approach in views matters rather differently in showing 
the equivalence of the SE to a pair of diffusion equations. The full theory is very 
elegant and extends to singular situations, etc. (cf. [5011^2). I* would be of interest 
here to further examine the quantum potential in this context. 

(5) In #6-#8 one arrives at a pair of equations by virtue of the "fractal" structure of 
space (where fractal here simply means that nondifferentiable paths are considered 
which generate a complex velocity). In |3j (as exhibited in ^,^8) one relates the 
quantum potential to the velocity showing its origin in the "fractal" derivative 
idea. 

We emphasize that in fact the quantum potential comes up in a serious manner in the Bohm 
theory, with refinements asin[niIl3II3l2n[l2lll2aiiSllllliSlllZlllHlEniE^ In 

fact, given that trajectories are at the base of this theory one can forsee a fractal Bohm 
theory in the future (cf . |671 1112j ) . On the other hand one can make convincing arguments 
for fields as the fundamental objects (except perhaps in the Bohmian type theories) with 
particules "emerging" (cf. |i64t 1142) ) as in quantum field theory (or perhaps via ripples or 
fractal structure in spacetime itself). 

It is not entirely clear how to handle derivatives in statistical or fractal theories. There are 
of course many powerful techniques available for Brownian motion and stochastic differential 
equations and there is a developing literature about differential calculus on fractals. Random 
walks and general discretization methods are also useful. Somehow one would like to imagine 
that the formal power of calculus (and duality via distribution like theories) might be strong 
enough to override the microscopic details about the domains of differential operators. 
Perhaps the coordinate derivative operators in situations such as #6-7^8 could be defined 
so that their domains are various fractal sets densely embedded in R" (in this connection 
see e.g. [201 El IZ3 EH Eni CSHl lESl 1113). in the end the most attractive formulation 
would seem to be some (more or less rigorous) version of a Feynmann path integral where 
precise definitions of the path space are not critical. 
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